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Abstract. Protostellar jets are present in the later stages of the stellar formation. Non-thermal radio 
emission has been detected from the jets and hot spots of some massive protostars, indicating the 
presence of relativistic electrons there. We are interested in exploring if these non-thermal particles 
can emit also at y-rays. In the present contribution we model the non-thermal emission produced in 
the jets associated with the massive protostar IRAS 18162-2048. We obtain that the y-ray emission 
produced in this source is detectable by the current facilities in the GeV domain. 
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INTRODUCTION 

Massive stars are formed within dense molecular clouds, accreting matter onto the 
central protostar with the formation of a circumstellar disc and forming bipolar jets 
and molecular outflows. In the last decade a handful of jets that emanate from massive 
protostars have been detected at radio wavelengths. In most of the cases, this radio 
emission has a negative spectral index (a), indicating the non-thermal nature of it. 
However, only in the case of the jet associated with the source IRAS 18162—2048, 
polarization of radio emission has been measured [1] confirming that the radiation is 
produced by the synchrotron process. 



The jets of the source IRAS 18162-2048 

The famous Herbig-Haro (HH) objects called HH80-81 are the south component of a 
system of several radio sources, located at a distance of ~ 1.7 kpc. The central source 
has been identified with the luminous (L* = 1.7 x 1O 4 L ) protostar IRAS 18162-2048, 
and HH80-N is the northern counterpart of HH80— 81. The velocity of the jet close to the 
protostar has been estimated to be as large as vi ~ 1400 km s" 1 [2], whereas HH 80 and 
HH 81 (located at z ~ 2.5 pc from IRAS 18162-2048) are moving with slower velocity, 
~ 150 km s , and HH80-N appears to be at rest. Radio observations [2] show that 
the central source has a spectral index a ~ 0.1, typical of free-free emission, whereas 
HH80— 81 and HH80-N are likely non-thermal sources, with a ~ —0.3. Between the 
central protostar and the HH-objects there is a chain of bright thermal and non-thermal 
knots. In particular, at a distance z p = 0.5 pc from the central protostar, an elongated 
structure of polarized radio emission has been recently detected through observations 



carried out with the Very Large Array (VLA) [1]. The detection of synchrotron radiation 
is evidence that there is a population of relativistic electrons in the source. These 
relativistic particles in the complex environment of the massive molecular cloud where 
the protostar is being formed can produce high-energy radiation through a variety of 
processes. In the present contribution we model the non-thermal emission produced in 
the polarized jet of IRAS 18162-2048, from radio to y-rays, at a distance z p = 0.5 pc 
from the central protostar, where the synchrotron emission is detected. 

Radio observations suggest that the jets of IRAS 18162—2048 are precessing, because 
the direction of them at z v is different than the direction where the HH objects are 
located. We assume that now the head of the jet is at z v , and HH80-81 and HH80 N 
are old ejections, ~ 4000 yr before. The head of the jet moving through the molecular 
cloud forms two shocks: a reverse shock in the jet, and a bow shock in the cloud, with 
velocities v rs |b s = vj — 3vb s /4 and vt> s = Vj/(1 + \J\j%), respectively, where % = n-]/n mc 
and Hi and n mc are the densities of the jet and the molecular cloud, respectively. At z p 
the width of the jet is 7?j ~ 10 17 cm, given nj = M/(7tRjvj) ~ 2.6 x 10 2 cm~ 3 , where 

M ~ 10~ 5 M yr" 1 [3]. At the location of HH80-N n mc (z ~ 2.5pc) ~ 400 cm~ 3 [4] 
is deduced from H a luminosity, and using that n mc ^z -1,5 [5] we obtain that at 0.5 pc 
« mc ~ 3 x 10 3 cm~ 3 . Thus we obtain that #(z p ) ~ 0.02 giving an adiabatic reverse shock 
with velocity v rs ~ vj and a radiative bow-shock. 

PARTICLE ACCELERATION AND NON-THERMAL EMISSION 

Being that adiabatic shocks are propitious to accelerate particles via the Fermi-I mecha- 
nism, we assume that the electrons (e) that emit the detected synchrotron radiation at z p 
are accelerated in the jet reverse shock (located also at z p ). Even though there is no ob- 
servational evidence of the presence of relativistic protons (p) in the jet of IRAS 18162- 
2048, these particles can be accelerated via Fermi-I as well as electrons, and we will take 
into account a population of non-thermal protons in our calculations. Particles acceler- 
ated in the reverse shock are injected in the downstream region following a power-law 
energy distribution Q e ^ p oc E~ 2 . The injection function Q is modified by radiative and es- 
cape losses in such a way that the part of the electron spectrum N e that produces the syn- 
chrotron emission between 1.5 and 15 GHz is N e E~ 2 . Considering equipartition be- 
tween the magnetic and the relativistic particle energy densities, i.e. B 2 / (8tt) = U e + U p , 
and the synchrotron fluxes given in [2], the magnetic field at z p results B ~ 0.1 mG [1]. 

Electrons lose energy via synchrotron, Inverse Compton (IC) scattering and relativis- 
tic Bremmstrahlung, but escape losses via downstream matter advection is more ef- 
ficient, with a timescale t^ y ~ 4R\/v\ ~ 4 x 10 9 (vj/(1000kms -1 )) -1 s as is shown 
in Figure 1 (left). Relativistic Bremsstrahlung losses have a timescale ?Brem ~ 1-4 x 
10 11 (nj/260cm~ 3 )~ 1 s, being this process the most efficient leptonic radiative channel 
at high energies. IC losses are not relevant because at z v the energy density of stellar pho- 
tons is low: £/ p h ~ L*/ (^z p c) ~ 9.4 x 10~ n erg s _1 . Relativistic protons radiate through 

proton-proton (pp) collisions on a timescale t pp ~ 1.9 x 10 11 (nj/260cm~ 3 )~ 1 s. Ad- 
vection losses constrain the maximum energy of electrons and protons given E e max = 
£p,max ~ 2.8 TeV, and N 6:P ~ Q e , p t adw . 
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FIGURE 1. Left: Acceleration, escape, and electron cooling timescales. The parameters of the jet 
and the molecular cloud used are specified in the plot. Right: Spectral energy distribution. Relativistic 
Bremsstrahlung and pp collisions are computed considering the shocked density of the jet (4nj = 1.1 x 
10 3 cm -3 ), and also the shocked density of the molecular cloud (n smc = 9x 10 s and 1.8 x 10 6 cm~ 3 , for 
vj = 700, and 1000 km s , respectively). 



Gamma-ray emission 

Knowing N e and N p we compute the spectral energy distribution considering the 
radiative processes mentioned in the previous Section (see Fig 1-right). We consider 
auto-absorption in the synchrotron spectrum and photon-photon attenuation at high 
energies, however this latter absorption process is not relevant at ~ 0.5 pc from the 
protostar. Considering two different values of vj: 700 and 1000 km s -1 , and the jet 
shocked density as target of relativistic Bremsstrahlung, and pp, the achieved bolometric 
luminosities of synchrotron radiation, IC scattering, relativistic Bremsstrahlung, and 
pp collisions result L sync ~5x 10 31 , L lc ~ 2 x 10 29 , LBrem ~ 8 x 10 30 , and L pp ~ 
5 x 10 31 erg s _1 , respectively. We note that in this case Vj only affect the values of 
E e ,max and E pmax , given very similar bolometric luminosities. 

We have mentioned that the most important cooling channels for electrons and protons 
are relativistic Bremsstrahlung and pp collisions, respectively. However, specific lumi- 
nosities not larger than ~ 10 30 erg s _1 (relativistic Bremsstrahlung) and 10 31 erg s _1 
(pp) are achieved in the GeV domain by these radiation mechanisms. In order to pro- 
duce a large amount of y-rays, an increment in the density of target particles is needed. 
This can occur if we consider that, being the large densitiy contrast between and mat- 
ter downstream the bow shock, both materials can be mixed by Rayleigh-Taylor (RT) 
instabilities. 

The radiative bow shock in the molecular cloud compresses the matter up to a density 
n smc ~ 10 5 (n mc /10 3 cm~ 3 )(vbs/ 100km s -1 ) 2 cm~ 3 after cooling via thermal emission 
up to reach a temperature of 10 4 K. Thus, being the large density contrast between 
4nj ~ 10 4 cm~ 3 and n smc , RT instabilities grow up to a wavelength l/i?j on a timescale 
t RT ~ 0.1?j, where fc ~ ZpAj [6,7]. Thus, at z p the emitter (i.e. the reverse-shock down- 
stream region) will be denser than 4rcj, and reaching a maximum value of n smc . Con- 



sidering the maximum increase in the density of the emitter, i.e. efficient RT mixing, 
detectable values of y-rays in the GeV domain are achieved, with specific relativistic 
Bremsstrahlung and pp luminosities larger than ~ 5 x 10 32 and 5 x 10 33 erg s _1 , respec- 
tively, as is shown in Fig 1 (right). However, note that detectable emission is achieved 
even in the case of less efficient mixing. By increasing the density of the emitter ~ 10 is 
enough to obtain detectable values of y-rays. 

CONCLUSIONS 

The massive protostar IRAS 18162-2048 presents very collimated and powerful jets, and 
an accretion disc, characteristics that put this source as the prototype of a massive star 
formation system. In addition to that, the jet of IRAS 18162-2048 has a large velocity 
(~ 1000 km s _1 ) and it is the unique proto stellar jet with detection of polarized radio 
emission, mapping the direction of the magnetic field parallel to the jet velocity. The 
large value of vj and the detection of polarization at radio wavelengths make the jets of 
IRAS 18162-2048 potential sources of y-rays [6]. 

In the present contribution we model the non-thermal emission (from radio to y- 
rays) produced at the location of the jet where synchrotron emission is clearly detected. 
Constraining the model with VLA data, and considering RT mixing between jet and 
molecular cloud shocked matter we obtain that the y-ray emission produced in the north 
jet of IRAS 18162-2048 can be detected by the Fermi satellite in the GeV domain. 
We note that besides IRAS 18162-2048, there are other massive protostars that can be 
also y-ray sources, as the case of IRAS 16547-4247 [8]. If high energy radiation from 
massive protostellar jets is detected, y-ray astronomy can be used to shed light on the 
star forming process, and cosmic ray acceleration inside molecular clouds. 
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